arXiv: 1503.05696v 1 [cs.IT] 19 Mar 2015 


Performance Analysis of Random Linear Network 
Coding in Two-Source Single-Relay Networks 

Amjad Saeed Khan and Ioannis Chatzigeorgiou 
School of Computing and Communications, Lancaster University, United Kingdom 
Email: {a.khan9, i.chatzigeorgiou} @lancaster.ac.uk 


Abstract —This paper considers the multiple-access relay chan¬ 
nel in a setting where two source nodes transmit packets to a 
destination node, both directly and via a relay node, over packet 
erasure channels. Intra-session network coding is used at the 
source nodes and inter-session network coding is employed at 
the relay node to combine the recovered source packets of both 
source nodes. In this work, we investigate the performance of 
the network-coded system in terms of the probability that the 
destination node will successfully recover the source packets of 
the two source nodes. We build our analysis on fundamental 
probability expressions for random matrices over finite fields and 
we derive upper bounds on the system performance for the case of 
systematic and non-systematic network coding. Simulation results 
show that the upper bounds are very tight and accurately predict 
the decoding probability at the destination node. Our analysis 
also exposes the clear benefits of systematic network coding at 
the source nodes compared to non-systematic transmission. 

Index Terms —Network coding, fountain coding, multiple ac¬ 
cess relay channel, decoding probability, block angular matrix. 

I. Introduction 

In their seminal paper Q~j, Ahlswede et al. proposed network 
coding as a means to share resources by allowing intermediate 
network nodes to combine packets prior to forwarding them 
to the next-hop nodes. Network coding, which can improve 
throughput in both wired and wireless networks, is classi¬ 
fied into intra-session and inter-session network coding El. 
The former approach mixes packets within the same flow, 
for example, packets of the same source node, and can be 
implemented using fountain codes. The latter scheme mixes 
packets of different flows, for example, packets that originate 
from different source nodes. 

Network coding can be combined with wireless node co¬ 
operation to reap the benefits of spatial diversity. A system 
configuration that is often considered in conjunction with 
network coding is the multiple-access relay channel (MARC). 
In the case of two source nodes, source packets are transmitted 
to both a relay node and a destination node on orthogonal or 
non-orthogonal channels. The relay node linearly combines the 
received packets of the two source nodes using inter-session 
network coding and forwards the network-coded packets to 
the destination node. The diversity gain and the performance 
of the system, measured in terms of the probability that the 
destination node will fail to recover the source packets of any 
of the source nodes, was studied by Bao and Li [3) and Chen 
et al. ffl. Woldegebreal and Karl 0 considered a similar 
setup to that in 0, 0, and investigated the impact of two 


different relay operating modes on the coverage area of the 
system. In an effort to improve the coding gain of the MARC, 
joint channel coding and network coding was introduced, for 
example in 0, 0, while compute-and-forward was proposed 
in ||0 as a means to increase the achievable rate of the MARC. 

An aspect that 0-0 have in common is that network 
coding is not used at the source nodes. Nevertheless, systems 
that employ intra-session network coding both at the source 
and the relay stages have been analyzed in the literature, 
for example |9l- lfTTil . However, the aforementioned work is 
concerned with the transmission of network-coded packets 
from a single source to a single destination via multiple relays. 
Consequently, the destination is not required to recover and 
identify source packets of different source nodes. 

This paper is concerned with transmission over an orthogo¬ 
nal MARC when two source nodes are present. The motivation 
of this work is to study a network configuration that encom¬ 
passes both intra-session network coding at the source nodes, 
as in ®-m and inter-session network coding at the relay 
node, as in 0-0. The probability that the destination node 
will successfully recover the source packets of both source 
nodes is used as the performance measure of the system. Our 
analysis builds on and extends recent work on random binary 
block angular matrices lfl2l . In our study, we have looked 
at the decode-and-forward relaying scheme, that is, network- 
coded packets received by the relay node are decoded and 
re-encoded before they are forwarded to the destination node. 
The derived probability expressions could be adapted to other 
network-coded relaying strategies that incorporate both intra¬ 
session and inter-session network coding schemes, as in fUl . 
or be used as benchmarks in performance comparisons. 

The rest of the paper is organized as follows. Section [IT] de¬ 
scribes the two-source single-relay system model in detail and 
formulates the problem statement. Section [III] builds on fun¬ 
damental probability expressions for decoding network-coded 
packets and obtains tight closed-form bounds on the proba¬ 
bility of the destination node recovering the source packets of 
both source nodes. Section[|V]validates the theoretical analysis 
and discusses simulation results. The contributions of the paper 
and future research directions are summarised in Section [V] 

II. System Model and Problem Statement 

We consider a network comprising two source nodes Si 
and S 2 , a relay node R and a destination node D, as shown 
in Fig. Q] Nodes Si and S 2 segment data into I\ 1 and K-> 


equally-sized packets, respectively. Let u \,..., uk 1 denote 
the source packets of node Si while UKi+i, ■ ■ ■ ,uk 1 +k 2 
represent the source packets of node S2. Each source node 
employs random linear network coding to combine source 
packets and generate coded packets. In non-systematic network 
coding, each source transmits Ne > Kp coded packets, where 
l = 1 , 2 . In systematic network coding, the first Kp transmitted 
packets are identical to the source packets, while the remaining 
Nf — Kf packets are coded. As is customary in network coding, 
each coded packet is transmitted along with a coding vector, 
which contains the Kp coefficients of the respective linear 
combination. In this paper, we consider coefficients that are 
chosen uniformly at random from the elements of the Galois 
field GF( 2 ). Therefore, each coded packet is the bitwise sum 
of source packets. 

Links between network nodes are modelled as packet era¬ 
sure channels. We use pe.T>, pe t R and pr.d to denote the 
packet erasure probabilities of the links connecting the f-th 
source node with the destination node, the Gth source node 
with the relay node and the relay node with the destination 
node, respectively. We assume that source nodes transmit on 
orthogonal channels enabling both the relay and the destination 
nodes to distinguish transmissions between the source nodes. 

The communication process is split into two phases. In 
the first phase, nodes Si and S2 transmit Ai and A2 coded 
packets, respectively, to node D. Node R overhears the trans¬ 
missions of the source nodes, stores the successfully received 
coded packets and attempts to decode them. Let Me and M' t 
be the number of coded packets from node Sf that were 
received by the destination node D and the relay node R, 
respectively. The coding vectors of the received coded packets 
can be stacked together at the receiving nodes to form coding 
matrices. At the end of the first phase, the coding matrices at 
nodes D and R can be expressed in block diagonal form as 
follows 


CsD 
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( 1 ) 



Fig. 1. Block diagram of a network consisting of two source nodes Si and 
S2. a relay node R and a destination node D. The packet erasure probability 
of each link as well as the number of transmitted and received coded packets 
at each node are also depicted. 


matrix at node D is 



which is a (Mi + M2 + Mr) x (AT 1 + K2) block angular 
matrix. Note that Crd has been expressed as the concatenation 
of matrices Cri and Cr 2, which were generated by node R 
and describe linear combinations of source packets originating 
from nodes Si and S2, respectively. We stress that the coding 
vector transmitted by node R is twice the size of the coding 
vectors transmitted by Si and S2; however, all coded packets 
in the network have the same size, which is customarily taken 
to be considerably larger than the size of the coding vectors. 

The objective of this paper is to characterise the system per¬ 
formance of the considered two-source relay-aided network. 
More specifically, we will carry out a performance analysis 
to determine the probability that the destination node D will 
decode the A'1 + K-> source packets of both nodes Si and 
S2, given that node D has recovered at least K\ + A'2 coded 
packets, that is, (Mi + M2 + Mr) > A'i + A'2. The impact of 
the chosen values for A and Ar on the system performance 
will also be discussed. 

III. System Performance Analysis 


where C/ is a Me x Kp matrix constructed at node D 
using the received coding vectors from node S^, and Cj is 
a Mp x Kp matrix that consists of the received coding vectors 
from node S/ : at node R. The dimensions of CgD and Csr 

are (Mi + M 2 ) x (Afi + A' 2 ) and (M{ + M') x (A'i + Af 2 ), 
respectively. 

In the second phase, if the relay node R successfully 
recovered the source packets of one or both source nodes, 
it linearly combines them in order to generate Ar coded 
packets. Therefore, the coding vector that accompanies each 
relay-generated coded packet consists of A'i + A'2 entries. 
If the relay node failed to decode the packets of either Si 
or S2 then the first A'i entries or the last A'2 entries of the 
coding vector, respectively, are set to zero. If Mr of the Ar 
transmitted coded packets are received by the destination node 
D, a Mr x (ATi + A'2) coding matrix Crd will be created and 
appended to CgD- At the end of the second phase, the coding 


Fundamental probabilities related to the rank of random 
matrices in GF( 2 ) are summarised in this section and are 
subsequently used in the derivation of expressions for the prob¬ 
ability that the destination node D will successfully recover the 
source packets of both source nodes, when they employ either 
non-systematic or systematic random linear network coding. 


A. Preliminaries: Fundamental Probability Expressions 


Let X be a m x k binary random matrix with m > k. 
We say that X is a full-rank matrix if the rank of X is k or, 
equivalently, k of the m rows of X are linearly independent. 
The probability of X being a full-rank matrix can be obtained 
as follows 


P (to, k) 


f( m i k) 

2 mk 


( 3 ) 


where 2 mfe is the number of all m x k binary matrices and 
f(m, k) is the number of all full-rank m x k binary matrices 






given by m 


k -1 

f(m,k)=] Q(2 m -2 l ). 

i=0 


( 4 ) 


The probability of X having rank r < k when m > r has also 
been derived in M and is equal to 


P r (m,k) = 2 ~ mk 


( f{m,r)f(k,r) \ 

V f{r,r) J 


( 5 ) 


Note that for r = k, expression 0 reduces to 10 - 

Let us now assume that matrix X has the following con¬ 
strained structure 


X = 



( 6 ) 


where the dimensions of submatrices A, B, C and D are 
a x a', b x b', c x a' and c x b', respectively. Matrices of 
this type, which are known as block angular matrices, were 
studied in El- It was proven that the probability of X being 
full-rank is given by 


P(o, a’, b,b',c )='%2 Pi(o, a')Pj(6,6')P(c, a’+b’-i-j). ( 7 ) 

i+j>a'-\-b' —c 

As implied by 0 , the rank of matrix X is a' + V if submatrix 
A has rank i, submatrix B has rank j and the remaining 
a' + b' — i — j columns of X are linearly independent, for all 
valid values of i and j. 

Expressions 0 . 0 and 0 will be invoked in the sub¬ 
sequent performance analysis. Note that character P is used 
exclusively to denote probabilities associated with the rank of 
matrices but character P is used to refer to probabilities related 
to the system model under consideration. 


B. Decoding Probability for Noil-systematic Network Coding 

In the general case of point-to-point communication over 
a channel with erasure probability p, the probability of the 
receiving node recovering all of the K source packets when N 
coded packets have been transmitted can be written as follows 

N 

P(N,K,p) = J 2 B{M,N,p)¥{M,K). (8) 

M=K 

The term B(M,N,p) denotes the probability mass function 
of the binomial distribution, that is, 

B(M, N,p) = -p) M p N - M . ( 9 ) 

Expression 0 enumerates all possible scenarios of retrieving 
the K source packets when M > K coded packets have been 
successfully received and have formed a full-rank M x K 
coding matrix. 

In the particular case of the considered relay-aided network, 
the probability that the destination node D will recover the 
source packets of both source nodes can be decomposed into 
the following three components: 


1 ) Unaided communication: Even though the relay node 
R has been deployed in the network, the destination node D 
could recover all of the source packets without the help of node 
R. The implies that both submatrices Ci and C2 in (0 are 
full-rank matrices and, consequently, Csd is also a full-rank 
matrix. Therefore, the probability that node D will recover 
the A'i + K2 source packets based solely on the N\ + Ng 
transmitted coded packets can be obtained using 0 as follows 

P s = P(N x ,Ki,p x , D ) P{N 2i K 2 ,P 2 ,u). ( 10 ) 

2 ) Partially aided communication: In this mode, the desti¬ 
nation node recovers the Kg source packets of node Sf based 
on coded packets transmitted both via the relay node and 
over the direct link between Sg and D. The destination node 
retrieves the source packets of the other source node, denoted 
by S^ where £ = 1,2 and £ 7^ £, without the assistance of 
the relay node. The probability that node D will recover the 
I\ 1 + A2 source packets, when transmission from node S/ : is 
aided by the relay node R while transmission from node 

is unaided, can be upper-bound by the following product 

Ps e HD < P{ N b K hPi,T>) P(N e ,K e , Pe , R ) 

N e 

■ V' B(Me,Ne,pe, d) 

M t = 0 ( 11 ) 

min(M^,/Q — 1) 

■ ^ Fi(Me, Ke)P(NR, Kt — *,Pr,d)- 
i=o 

The first two terms on the right-hand side of (fill represent 
the probability that nodes D and R will recover the source 
packets of nodes S/ and S*, respectively, when the direct links 
are used. The remaining two lines compute the probability that 
node D will construct a coding matrix of rank Kg by obtaining 
i linearly independent coding vectors from node Sf and Kg — i 
linearly independent coding vectors from node R. Derivation 
of this probability invoked and extended a degraded version 
of the right-hand side of 0 , where X in 0 was redefined as 
X = (A C) T . 

The reason that the right-hand side of ifTTl) is an upper bound 
and not the exact expression for Ps,,rd lies to the fact that 
the probability of the relay node decoding the packets of node 
S/ : is not independent of the probability that the destination 
node will decode the packets of the same node. For example, 
consider the case when Ng = 10 coded packets are transmitted 
to both D and R and pg,\g = pgj{ = 0 . 1 . If each node recovers 
9 coded packets then each node will have at least 8 of them 
in common. Therefore, if node D fails to recover the source 
packets of node Sf, node R will most likely also fail to recover 
them and will not be in the position to assist node Sf in its 
transmission. However, as the value of the product Ng pg jg or 
Nipg ,d increases, the upper bound gets tighter, as will become 
evident in Section [IV] 

Using tlTTb . the probability that the destination node will 
recover the source packets of both Si and S2, when either Si 
or S2 is aided by the relay node R, is given by 

-PsRD = PsiRD + Ps 2 RD- 


( 12 ) 



3) Fully aided communication: In this case, both Si and 
S 2 need the aid of the relay node R in order to deliver the 
necessary number of coded packets to the destination node. 
Node D successfully decodes the coded packets transmitted 
via node R and over the two direct links, and recovers all 
source packets. The probability that node D will recover the 
K\ + K ‘2 source packets, when both source nodes are assisted 
by the relay node, can be upper-bound as follows 

-Prd < P{Ni, P(N 2 , K 2 ,p 2 , R ) 

tVi n 2 

■ ^ B(Mi, 7Vi,pi,d) ^ B(M- 2 , N 2 ,p 2 ,n) 

Mi=0 M 2 = 0 

l max Jmax 

•EE 

i =0 j —0 

• P(N r, Ki+K 2 —i—j, pr,d)- 

The first line on the right-hand side of (fl3l > expresses the prob¬ 
ability that node R will recover the source packets of both Si 
and S 2 . The remaining lines compute the probability that node 
D will receive i, j and K\ + K 2 — i — j linearly independent 
coding vectors from Si, S 2 and R, respectively, for all valid 
values of i and j. Similarly to (fill , we set the upper limit of 
the third sum in ( IT3l > equal to i max = min(Mi, K 1 — 1); this 
ensures that the number of linearly independent coded vectors 
i, which have been received directly from node Si, is neither 
greater than the total number of received coded vectors Mi, 
nor equal to or greater than the number of source packets I\\. 
The definition of ?' rnax prevents i from taking the value I\\ 
because cases where node D can recover the K\ source pack¬ 
ets without the help of node R have already been considered 
in unaided and partially aided communication. Following a 
similar line of reasoning, we set j max = min(M 2 , K 2 — 1) 
in (0. Observe that the last two lines of ( fl3l > constitute a 
formula that is a constrained extension of (|7). 

The overall decoding probability at the destination node D 
can be obtained by adding the three constituent probabilities, 
that is, 

Pd = Ps + Psrd + Prd- (14) 

We remark that if the right-hand side of ( | I Ok (I I 1 1 and ( fTTl i 
are used in ( 1 1 4t i to compute Pg, Psrd and Prd, respectively, 
an upper bound on Pd will be obtained. 


C. Decoding Probability for Systematic Netw’ork Coding 

In H5), systematic network coding for point-to-point com¬ 
munication was studied and it was proven that the probability 
of a receiving node decoding all of the K source packets, given 
that K < M < N packets have been successfully received, is 


P'{M,K, N) 


( K \ ( N ~ K \ 
2-jh=h mi n\h)\M-h> 


F(M-h,K-h ) 



(15) 


where h m i n = max(0, M — N+K). Expression ( IT5l > considers 
the possibility of receiving h systematic and, hence, linearly 
independent packets out of the K transmitted systematic pack¬ 
ets and computes the probability that there exist K — h linearly 


independent coded packets among the remaining M — h 
packets, for all valid values of h. Following the same line 
of reasoning as in Da, we can express the probability of 
receiving r < K linearly independent coded packets as 


Pf(M,K,N) = 


EL b^QC^Pr-biM-h^K-h) 

(m) 


(16) 


provided that M > r. Similarly to the case of non-systematic 
network coding, the probability of the receiving node recov¬ 
ering all of the K source packets when N packets have been 
transmitted, denoted by P'(N, K. p), can be obtained from ([8]) 
by replacing P (M,K) with P'(M, K, N). 

Taking into account (IT5l) and (IT6l) and using the same train 
of thought as in Section IIII-BI we can obtain an expression 
for the performance of the considered two-source single-relay 
network for the case of systematic network coding. More 
specifically, the probability that the destination node will 
recover the source packets of both source nodes is given by 


Pn = Ps 


(Ps 


SiRD 


Ps 2 rd) + -Prd 


(17) 


where 

Ps = P'iN^K^p^) P'(N 2 ,K 2 ,p 2 v), 

Ps.RD < P’(N h K h p lD ) P’(Nt,K eM ) 

N e 

■ ^ B(Me,Ne,pep) 

Mg =0 

min(Mg,Kg — 1) 

■ ^ Pl(M e ,K e ,N e )P(N R ,K e -i,p R>n ) 

i =0 

for £ = 1,2, and 

P R D < P'(iVi, K\,p rr ) P'(N 2 ,K 2 ,p 2tR ) 

N± N 2 

■ ^ B{Mi,Ni 7 pijf) ^ B(M 2 ,N 2 ,p 2 ,v) 

Mi=0 M 2 = 0 

Imax Jmax 

■EE Pl(M u K u NJPftMi, K 2 , N 2 ) 

i —0 j =0 


( 18 ) 


(19) 


( 20 ) 


• P(N r , Ki+K 2 — i—j , pr,d)- 

The validity and tightness of the derived performance 
bounds will be investigated in the following section. 


IV. Results and Discussion 
In this section, comparisons between the derived theoreti¬ 
cal upper bounds and simulation results will be carried out 
for both systematic and non-systematic network coding. For 
convenience, a symmetric network configuration has been con¬ 
sidered, according to which K\ = K 2 = K, N\ = N 2 = N, 

Pi,d = P2,d = Ps,d and pi, R = p 2 ,R = Ps,R- 

Fig. [2] compares simulation results with the theoretical 
expression in (IT4l) as a function of N R , for different values of 
K and N. As explained in Section llll-BL the interdependency 
between the decoding probability at node R and the decoding 
probability at node D is evident when K = 10 and N = 15; in 
this case, the upper bound yields a marginally higher decoding 









Fig. 2. Comparison between theoretical upper bounds obtained from GD and 
simulation results for different values of K and N. The erasure probabilities 
have been set to pg D = 0.3, Ps R — 0.1 and pr d = 0.2. 

probability than that obtained via simulations. However, the 
interdependency becomes smaller and the upper bound gets 
tighter with an increasing number of source packets K and, 
consequently, an increasing number of transmitted packets N. 
We observe that for K = 20 and N = 30, the derived upper 
bound coincides with the simulation results. 

The tightness of the proposed upper bound is also illustrated 
in Fig. [3] which depicts the impact of the source-to-destination 
channel quality, represented by ps.n, and the number of coded 
packets Nr transmitted by the relay node on the system 
decoding probability Pd- As expected, aid by the relay is of 
key importance to the source nodes as the quality of the direct 
channel between each source node and the destination node 
deteriorates. The theoretical bounds accurately quantify the 
relationship between pg.D and the number of coded packets 
Nr that need to be transmitted by the relay to achieve a target 
decoding probability. 

Fig. E carries out a performance comparison between sys¬ 
tematic and non-systematic network coding (NC) for various 
values of ps,R- As is evident from the figure, if systematic 
NC is used at the source nodes and the source-to-relay 
channel conditions are good, the destination node requires 
fewer excess coded packets N — K from the source nodes 
to correctly decode all of the K\ + K2 source packets. This 
observation is in agreement with the findings in lfT5l for 
point-to-point communication. As the source-to-relay channel 
quality deteriorates, systematic NC performs similarly to non- 
systematic NC. Nevertheless, systematic NC still offers the 
benefits of progressive packet recovery and reduced decoding 
complexity, as detailed in 02). 

V. Conclusions and Future Work 

This paper studied the performance of a network comprising 
two source nodes transmitting to a destination node via a relay 
node, where random linear network coding is used both at 
the source nodes and the relay node. Upper bounds on the 
probability of the destination node successfully recovering the 
packets of both source nodes were derived for both systematic 
and non-systematic network coding. Simulation results con¬ 
firmed the validity of our theoretical analysis and established 
that the upper bounds get tighter and accurately predict the 



Fig. 3. Comparison between theoretical upper bounds obtained from 
and simulation results for different values of ps,r >• The remaining system 
parameters have been set to K = 20, N = 30, Ps R — 0-3 and pr,d — 0.2. 



Fig. 4. Performance comparison of systematic and non-systematic network 
coding as a function of the excess coded packets N — K transmitted by each 
source node for various values of ps,R- The remaining system parameters 
have been set to K = 20 , N r = 15 , Ps D = 0-3 and pr d — 0 . 1 . 

system decoding probability for an increasing number of 
transmitted coded packets by the source nodes. Furthermore, 
we demonstrated that systematic network coding can yield 
a similar or better performance than non-systematic network 
coding depending on the quality of the uplink channels. 

Future work will extend the system model to networks of 
multiple source nodes and, possibly, multiple relay nodes, 
and will aim to obtain expressions for the overall decoding 
probability. In this endeavour, both binary and non-binary 
network coding will be considered, that is, coded symbols in 
transmitted packets will be elements of a finite field GF(g) 
of size q > 2. An additional line of investigation to be 
pursued in the short term is the optimisation of the ratio 
between the coded packets transmitted by each source node 
and the coded packets transmitted by the relay node in order 
to achieve a target decoding probability, while meeting the 
energy constraints of the source nodes. We will also strive 
to draw comparisons of the considered decode-and-forward 
scheme with other network-coded forwarding protocols as well 
as with networks employing physical network coding. 
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